Thermal phase transitions in the attractive extended Bose-Hubbard Model with 
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By means of quantum Monte Carlo simulations implemented with a two-loop update scheme, 
the finite-temperature phase diagram of a three-body constrained attractive Bose lattice gas is 
investigated. The nature of the thermal phase transitions around the dimer superfiuid and the 
atomic superfiuid is unveiled. We find that the Z2 symmetry-breaking transitions between these 
two superfiuid phases are of first order even at nonzero temperatures. More interestingly, the 
thermal transition from the dimer superfiuid to the normal fluid is found to be consistent with the 
Kosterlitz-Thouless type but giving an anomalous universal stiffness jump. It demonstrates that 
this transition is driven by unbinding of pairs of fractional vortices. 
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Since the early proposals of possible pair superfluid- 
ity in an attractive Bose gas?i the existence of this state 
of matter at low temperatures has been intensively ex- 
plored. Recent developments on the manipulation of 
ultracold gases may provide direct routes in the search 
for the pair superfiuid phase and in probing quantum 
critical behaviors around this intriguing state. In the 
context of Bose gases in the continuum, the possibility of 
observing such a pair supcrfluid phase near a Feshbach 
resonance has been proposed4r"— Another realization in 
attractive bosonic lattice gases with three-body on-site 
constraint is also suggested recently^ While the discus- 
sions of this pair condensed state was originally focused 
on the system of bosonic particles, similar physics can be 
applied to explain some exotic phases in other physical 
systems. For example, the quantum spin nematic state 
in some frustrated spin systems can be understood as the 
condensate state of bound magnon pairsi^^i^ 

The pair superfiuid phase or the dimer superfiuid 
(DSF) phase consists in the formation of a macroscopic 
coherent state made of boson pairs. Here we focus on the 
case of single-species bosons. Contrast to the conven- 
tional atomic superfiuid (ASF) with non-vanishing or- 
der parameters (a) and (a^) 7^ (here a denotes 
the boson annihilation operator), DSF is characterized 
by the vanishing of atomic order parameter ((a) — 0) 
but nonzero pairing correlation ((a^) ^ 0). Apart from 
the above local order parameters, one can use superfiuid 
stiffness to identify the superfiuid states. It is expected 
that, in the DSF phase, pairs of bosons will wind to- 
gether around the system, such that only even numbers 
of winding can occur. Hence the DSF phase can also 
be characterized by a nonzero pair superfiuid stiffness 
given by even winding numbers pcvon 7^ and a zero 
stiffness corresponding to odd winding numbers podd = 
[see Eq. (I2])]iii Such an even-odd effect will not appear 
in the ASF phase. From the symmetry-breaking perspec- 
tive, the U{1) symmetry of the global phase transforma- 
tion a — ;> e^'^'a is completely broken in the ASF phase, 
while a residual Z2 discrete symmetry remains in the DSF 
phase due to the 7r-periodicity of (a^). As a result, there 
should exist an Ising-like quantum phase transition be- 



tween these two phases upon tuning system parameters. 
Nevertheless, as pointed out in Refs. MM, quantum fiuc- 
tuations can turn this transition into a first-order one due 
to the Coleman- Weinberg mechanismji^ Although there 
are many theoretical investigations on the existence of 
the DSF phase at zero temperature and on the nature 
of the related quantum phase transitions^"— i"'^^'^^ reli- 
able quantitative predictions, especially for the DSF- ASF 
transition, have not yet been provided except for the one- 
dimensional caseJ^ Moreover, the physics of the thermal 
transitions out of the DSF phase has neither been ad- 
dressed. 

In the present work, the nature of the finite- 
temperature phase transitions around the DSF phase in 
attractive bosonic lattice gases with three-body on-site 
constraint^ is explored numerically. We employ here the 
stochastic series expansion (SSE) Monte Carlo methodic 
generalized by allowing pair updates such that two in- 
dependent loops can merge and move together. The im- 
portance of two-loop (or two-worm) algorithms for effi- 
cient sampling of the DSF phase (or any paired phase) 
in large system sizes has been discussed in the litera- 
turcjiiii^"— Our results for the phase diagram are sum- 
marized in Fig. [TJ The existence of the DSF phase char- 
acterized by the even-odd effect of the stiffness is con- 
firmed in our simulations. Besides, it is found that the Z2 
symmetry-breaking transitions between the DSF and the 
ASF phases are still of first order at finite temperatures. 
This indicates that the underlying Coleman- Weinberg 
mechanism for the DSF- ASF transition at zero temper- 
ature is not completely spoiled by the thermal fiuctua- 
tions. As for the transitions to normal fiuids (N), both 
the DSF-N and the ASF-N transitions are found to be 
of continuous Kosterlitz-Thouless (KT) type<^ Remark- 
ably, these two KT transitions have distinct characters. 
Our data support an anomalous value of universal stiff- 
ness jump at the KT transition out of the DSF phase, 
which is four times larger than that out of the ASF phase. 
This observation clearly establishes that the DSF-N tran- 
sition is driven by the proliferation of unconventional 
topological defects (half- vortices). Thus the anomalous 
stiffness jump can serves as a unique signature for the 
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FIG. 1: Finite-temperature phase diagram of a three-body 
constrained attractive Bose gas on a square lattice as de- 
scribed in Eq. Q with |;7| = 1, 1/ = 0.25, and yL = -0.55. 
The solid lines indicate the DSF-N and ASF-N transitions of 
continuous KT type, and the dashed line shows the first-order 
DSF-ASF transition. Lines are guide to eyes. 

search of the DSF phase in experiments. 

We consider the following extended Bose-Hubbard 

t 3 

model with a three-body constraint a\ =0, 

+ V^^ninj - fi^^m . (1) 

Here, ai{al) is the bosonic annihilation (creation) oper- 
ator at site i, t is the hopping matrix element, U < 
the on-site two-body attraction, and ^ the chemical po- 
tential. V > denotes the nearest-neighbor repulsion, 
which can come from the dipole-dipole interactions of 
the dipolar bosons polarized perpendicularly to the lat- 
tice plane by truncating it off at the nearest-neighbor dis- 
tance. The convention signifies a sum over nearest- 
neighbor sites i and j. The on-site constraint can arise 
naturally due to large three-body loss processes?^ and it 
stabilizes the attractive bosonic system against collapse. 
The intriguing quantum critical behaviors for y = at 
zero temperature have been explored in Refs. [sifisl Here 
we consider a finite repulsive V which is found (not shown 
here) to stabilize the DSF phase and extend the region of 
this phase to larger t. It thus facilitates the access of the 
DSF phase in our simulations, and may also allow easy 
observation in the real experiments as well. We have also 
confirmed the same behaviors of the DSF for V = 0. As 
an illustration, we choose the parameters V = 0.25 and 
/i = —0.55 {\U\ = 1 as the energy unit), so that the boson 
densities are always smaller than half-filling. 

To characterize different phases, several kinds of super- 
fluid stiffness are evaluated. In SSE, the conventional su- 
perfluid density ps at temperature T is computed by mea- 
suring the fluctuation of the winding number W within 
the simulations r^iS^ ps ~ mT{W'^). Here m = l/2t is 
the effective mass of the bosons in a square lattice. If 



a macroscopic fraction of the bosons winds around the 
system, the system will give a finite ps- In the ASF 
phase, the usual algorithm is able to let a large number 
of bosons wind around the system and is known to be 
efficient. However, in the DSF phase, bosons are paired 
into dimers and hop together as pairs. Therefore only 
even winding numbers can occur. If we define the super- 
fluid stiffness Pcvon (odd) with respect to the even (odd) 
winding number Wcvon (odd) in the following way: 

Pcvcn (odd) = "l?'((W^GVcn (odd)) ) , (2) 

then a clear even-odd effect, where ps — Povon ^ and 
Podd = 0, should be observed in the DSF phasCf^i It is 
worth to note that there is no pair hopping term in the 
Hamiltonian H, so that the motion of boson pairs is man- 
ifested as a second-order effect in the single-particle hop- 
ping parameter t. Within the conventional one-loop up- 
dates of the SSE simulations, the boson numbers can only 
be varied by one during each vertex update. Here in our 
two-loop algorithm, the pair updates that change the bo- 
son number by two are included. As a consequence, two 
independent loops can merge and move together in such 
a way that effectively simulates the pair hopping in the 
DSF phase. Because the accumulation of a large num- 
ber of winding boson pairs is exponentially suppressed 
by the short-range nature of the usual one-loop updates, 
the two-loop update scheme is necessary to improve the 
efficiency of the algorithm in finding the DSF phase. 

As pointed out in Ref. @, when the attraction U is 
strong enough, pairing correlation among bosons can be 
nonzero such that the system locates in a DSF phase. 
Conversely, by lowering the ratio of \U\/t, the ASF state 
can be stabilized. Thus there is a transition separat- 
ing the DSF and the ASF phases. At zero temperature, 
this Z2 symmetry-breaking transition is shown to be of 
first order in most conditions.— By performing quantum 
Monte Carlo simulations to study the model in Eq. ^ 
on square lattices of size = L x L under periodic 
boundary conditions, we find that this conclusion is still 
true at low temperatures. The results of various kinds 
of superfluid stiffness defined above as functions of the 
hopping parameter t at temperature T = 0.005 with 
system size L — A8 are shown in Fig. [2l It is found 
that, Ps — Pcven 7^ and podd = for smaller hop- 
ping parameters, whereas for larger hopping parameters 
Ps ~ 2podd 7^ and there is no even-odd effect. This 
indicates a transition from the DSF with Podd = to the 
ASF with nonzero Podd- As seen from the lower panels 
of Fig. [21 upon increasing system sizes, podd shows an 
abrupt jump at i = — 0.121. It implies that the DSF- 
ASF transition remains being of first order at finite tem- 
peratures. To further exclude the possibility of a second- 
order transition, we have performed a finite-size scaling 
by assuming the critical exponents of the Ising univer- 
sality class. The data of different system sizes do fail to 
collapse into a universal curve. This again advocates that 
the DSF-ASF transition should be of first order. Same 
conclusion about the nature of the DSF-ASF transition 
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FIG. 2: (Color online) Upper panel: several kinds of super- 
fluid stiffness as functions of f at T = 0.005 with system size 
L = 48. Lower panels: data of podd for different system sizes 
are shown to demonstrate the finite size effects. All error bars 
are smaller than the symbol size if not shown. 



is reached for lower temperature T = 0.001. 

As seen in Fig. [TJ increasing temperature further, 
either the DSF or tlie ASF order will be eventually 
destroyed by thermal fluctuations, and will undergo 
a symmetry-restoration transition to the normal-fluid 
state. Before presenting the detailed analysis of these 
transitions, let's begin with discussions on some general 
aspects of their nature. While both the DSF-N and the 
ASF-N transitions are expected to be the continuous KT 
transitions, there is an essential difference between them. 
Because the DSF phase preserves the tt phase-rotation 
symmetry, it can be characterized by an algebraic order 
in cxp(2i0) rather than in cxp(i0). where 9 is the super- 
fluid phase. Therefore, the KT transition out of the DSF 
phase should consist in proliferating of pairs of fractional 
vortices with vorticity v = ±1/2, instead of the ordinary 
(integer) vortices. As a result, the strength of the loga- 
rithmic interaction of these fractional vortices will be re- 
duced by a factor oiv^ = 1/4 in comparison with that of 
integer vortices, and the KT transition temperature Tkt 
will be decreased substantially. The universal jump of the 
superfluid stiffness at such KT transitions driven by un- 
binding of fractional vortices has an anomalous value 
Ps = {2 / Trv'^)mTKT , which amounts to ps = SuiTkt/'!^ 
for vorticity v = ±1/2. The existence of an anomalous 
KT transition has been proposed in other physical sys- 
tems i^i^ 

It is known that the KT transition temperatures can 
be determined with good accuracy by utilizing the renor- 
malization flow and the universal jump of the superfluid 
density at the transition pointi ^^'^^ Here we generalize 
the data analysis suggested in Ref. [2^ to take into account 
the anomalous superfluid stiffness jump. Consistent re- 
sults can be reached by other kinds of analysis. Define 
R = nv^ ps/2mT such that i? = 1 at the KT transition 
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FIG. 3: (Color online) Superfluid density ps vs temperature T 
with different system sizes for t = 0.14. The dash line denotes 
the location of Tkt- Inset: solutions of the Eq. Q for different 
pairs of system sizes (Li, L2) for R = -kv^ pa/2mT with v — 1. 
The blue solid line is the hnear fit k = 1 -I- 8. 71 (Tkt - T)/t 
with Tkt — 0.038 and the red dash line is k = 1. 

temperature Tkt. It is known that the KT renormaliza- 
tion group equations can be cast into an integral form, 

41n(L2/ii) = r f , (3) 

Jr2 t^iMt) - i^) + t 

where the parameter n is an analytic function of temper- 
ature, and the temperature giving k, ~ 1 corresponds to 
the KT transition point. For T < Tkt, a- hnear function 
in (Tkt - T) is expected, that is, k(T) ~ l-f c(Tkt - T) 
with a positive slope c— Thus, after taking different pairs 
of system sizes in Eq. ([3]) at each temperature to deter- 
mine the k{T) curve, the location of the KT transition 
temperature can be achieved by finding k(Tkt) = 1- 

Our findings of the superfluid density ps as a function 
of temperature T for t = 0.14 are depicted in Fig. |31 
We can sec the strong system size dependence character- 
istic to the KT transition especially around and above 
the critical temperature. As seen from Figs. [1] and [21 the 
low-temperature states at this hopping parameter belong 
to the ASF phase. Thus the conventional KT transition 
driven by unbinding of ordinary (integer) vortices is ex- 
pected. The results of the parameter k in Eq. ([3]) with 
1/ ~ 1 ioi R are shown in the inset of Fig. [3] The val- 
ues of k{T) extracted from various pairs of system sizes 
clearly collapse into a straight line around the value of 
K = I. This smooth analytic behavior of k{T) supports 
that the transition is of the usual KT type. Besides, the 
KT transition temperature is found to be Tkt — 0.038. 
Thus the corresponding universal stiffness jump in the 
thermodynamic limit is ps = 2mT}<iT/'K ~ 0.086. 

As discussed above, the low-temperature states for 
small hopping parameter t can be in the DSF phase, and 
the anomalous KT transition driven by unbinding of half- 
vortices ~ ±1/2) should be observed when tempera- 
ture is increased. As an illustration, our results of ps vs 
T for t = 0.1, where the system stays in the DSF phase 
at low temperatures, are shown in the left panel of Fig. 31 
Similar to the previous data analysis but using v = 1/2 
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FIG. 4: (Color online) Left panel: superfluidity stiffness ps 
vs T with different system sizes for t = 0.1. The dash line 
denotes the location of Tkt- Right panel: solutions of the 
Eq. ([3} for different pairs of system sizes (Li, L2) for R — 
ps /2mT with v = 1/2. The blue solid line is the linear 
fit K = 1 + 10.67(rKT - T)/t with Tkt ^ 0.0036 and the red 
dash line is «: = 1. The inset shows the same analysis but 
using u = 1 for R. 



for i?, the results of k{T) are presented in the right panel 
of Fig. m Again, the data collapse is evident, support- 
ing that this transition is of the KT type, but originated 
from the unbinding of half- vortices. For comparison, we 
repeat the same analysis with v = 1 for R, as presented 
in the inset of the right panel of Fig. 01 Data collapse 
becomes poor around the value of k = 1, and n{T) docs 
not behave as a linear function in (Tkt — T) for T < Tkt- 
These observations eliminate the possibility of the con- 
ventional KT transition. Therefore, our findings offer 
strong evidences of the existence of a half- vortex unbind- 
ing transition out of the DSF phase and also serve as a 



clear signature of the DSF phase itself. From Fig. |4l we 
observe that the KT transition happens at Tkt — 0.0036, 
which is about one order less than the value for the case 
of t = 0.14. The corresponding stiffness jump for t = 0.1 
in the thermodynamic limit is ps = SmTKT/'"' — 0.046, 
which is again less than the value for t ~ 0.14. 

In our simulations, the above conclusions about the 
nature of the DSF-N and the ASF-N transitions apply 
also to the corresponding cases with the chosen i's as 
presented in Fig. [TJ This implies that, if it exists, the 
possible region of the superfluid-to-normal-fluid transi- 
tions of types different from the KT transition should be 
very narrow. 

In summary, we study the DSF and the ASF phases 
in the extended Bose-Hubbard model with a three-body 
constraint and identify the nature of the related thermal 
phase transitions. The even-odd effect in the stiffness dis- 
tinguishes the DSF from the ASF phase. The transition 
between these two superfluid phases is found to be of first 
order. More importantly, the DSF-N transition is shown 
to be driven by the topological defects of half-vortices, 
in great contrast to the conventional ASF-N transition. 
Stimulated by the recent progress on detecting the KT 
transitions in cold atom, — our predictions may be veri- 
fied in the near furture. 
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